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^ ; Abstract 

-(— > , 

^ ■ Rayleigh waves are considered for crystals possessing at least one 

plane of symmetry. The secular equation is established explicitly for 
'^ I surface waves propagating in any direction of the plane of symmetry, 

^ ' using two different methods. This equation is a quartic for the squared 

Q ■ wave speed in general, and a biquadratic for certain directions in cer- 

tain crystals, where it may itself be solved explicitly. Examples of such 
materials and directions are found in the case of monoclinic crystals 
^ ■ with the plane of symmetry at X3 = 0. The cases of orthorhombic 

^«. , materials and of incompressible materials are also treated. 
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1 Introduction 

The simplest physical setting involving a boundary value problem for an elas- 
tic material is that of a semi-infinite body with a plane boundary left free 
of tractions. However, the consideration of small amplitude deformations 
and motions of such a half-space leads usually to considerable mathematical 
difficulties, especially when the material is anisotropic. Indeed in the case 
of general (triclinic) anisotropy, the equations of motion (or of equilibrium) 
lead to the resolution of a sextic for the partial inhomogeneous plane waves 
(or deformations), whose roots cannot be obtained explicitly. Consequently, 
closed-form solutions have been sought for materials with at least orthorhom- 
bic symmetry, because then the equations of motion lead to a biquadratic, 
and because this case covers 16 different types of common symmetry classes 
such as tetragonal, hexagonal, or cubic (Royer and Dieulesaint, 1984). 

In between the classes of triclinic crystals (no plane of symmetry) and 
of orthorhombic crystals (three orthogonal planes of symmetry), is the class 
of monoclinic crystals, with only one plane of symmetry. Among the three 
possibilities for the orientation of the symmetry plane, the configuration of 
a half-space X2 > made of monoclinic material with the plane of symmetry 
at 0:3 = is particularly important for two-dimensional deformations, (a) 
because in-plane stress and in-plane strain decouple from anti-plane stress 
and anti-plane strain, respectively, so that the equations of motion yield a 
quartic; and (b) because these materials are structurally invariant (Ting, 
2000) that is, the stress-strain relationships retain their form with respect to 
rotations in the (xi, a;2)-plane about the xs-direction, so that results obtained 
along the material axes Xi, X2, x^, are easily transposed along the rotated 
axes xl, X2, X3, say. On the other hand, the problem of surface Rayleigh 
waves is also of prime importance because it is relevant to the study of many 
other problems for anisotropic elastic half-spaces, such as: near-the-surface 
stability analysis of a deformed half-space (Biot, 1965), normal forces applied 
to a half-space (Lamb, 1904), punch and indentation of half-space (Green and 
Zerna, 1968), steady state crack propagation (Broberg, 1999), and so on. 
An up-to-date account on the research and applications of surface acoustic 
waves in materials science can be found in (Hess, 2002). In the present 
paper, the secular equation is established for surface waves in monoclinic 
crystals with the plane of symmetry at x^ = 0, where by 'secular equation' 
is meant the function of the squared wave speed c^ which is zero when the 
tractions on the plane 0:2 = and at a;2 — )■ 00 are zero. This equation is 
valid for the propagation of a Rayleigh wave in any direction of a symmetry 
plane for crystals possessing one plane of symmetry and of course for crystals 
of a higher order symmetry class, such as orthorhombic symmetry. It is 



also easy to take an eventual incompressibility of the elastic half-space into 
account (Destrade et al., 2002). Finally, after a rotation about the a:3-axis, 
the quartic secular equation may reduce to a biquadratic which can then be 
solved explicitly. 

The secular equation is obtained in Section 3, after the equations of mo- 
tion and boundary conditions for the problem have been recalled in Section 
2. This equation is obtained in two different manners, first in "covariant" 
(Furs, 1997) form, then as a quartic in the squared wave speed (Currie 1979, 
Destrade 2001, Ting 2002a). In the final Section (§4), the results are applied 
to other situations. First, a rotation is made for the (xi,X2) plane about 
the Xa-axis and, at least for three monoclinic crystals (diallage, gypsum, tin 
fluoride), two directions are found for which the secular equation may be 
solved explicitly. Then the results are specialized from monoclinic to or- 
thorhombic symmetry. Finally, the constraint of incompressibility is taken 
into account, and a numerical problem left open by Nair and Sotiropoulos 
(1999) is resolved. 

Throughout the paper, the dynamical analysis is based on the use of the 
components of the tractions rather than the displacements, and expressions 
are found in terms of the stiffnesses as well as in terms of the reduced com- 
pliances. 

2 Preliminaries 

Here we recall the equations of motion for a linearly elastic semi-infinite 
body, made of a monoclinic material with the plane of symmetry at 2:3 = 0, 
and seek a solution in the form of a surface wave solution, that is a solution 
which propagates in the Xi-direction, leaves the plane X2 = free of tractions, 
and vanishes as X2 — ?■ c>o. Because for such materials, in-plane motions are 
decoupled from anti-plane motions (Stroh, 1962), it is sufficient to seek a 
solution in the form of a two-component displacement vector u, such as 

u(xi,X2,t) = [f/l(x2),t/2(x2),0]V^(^-^*\ (1) 

where f/i, and f/2 are functions of a;2 satisfying f/i(oo) = f/2(oo) = 0, A; is the 
wave number, and c is the wave speed. 

With this convention, the equations of motion are written as (Mozhaev, 
1995) 

aU" + i/SV - 7U = 0, (2) 

where U = [f/i,f/2]^ and the prime denotes differentiation with respect to 
kx2- Here the symmetric 2x2 matrices aij, Pij, and 'jij, are given in terms 



of the elastic stiffnesses C's and of the mass density p by 
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Finally for the problem at hand, the following boundary conditions must also 
be satisfied: 



Ce,U[iO) + C26f/^(0) + ^Ci6f/i(0) + tC,eU2iO) = 0, 

C26f/((0) + C22U!,iO) + tCuUliO) + 2C26f/2(0) = 0. 



(4) 



Dual to this approach is one involving the components of the tractions 
acting upon the planes parallel to the free surface, instead of the components 
of the mechanical displacement. Indeed, just as in-plane strain is decoupled 
from anti-plane strain, so is in-plane stress from anti-plane stress (Stroh 1962, 
Ting 1996, Destrade 2001a). Thus, introducing the scalars functions ti{x2) 
and ^2(3^2), defined by 

cr2l{x^,X2,t) = h{x2)e'''^'-''-''\ a22{xi,X2,t) = t2(x2)e^^(^^-^*\ (5) 

where o"2i and (722 are the in-plane stress components, the equations of mo- 
tions may be written as (Destrade, 2001) 



at" - 0i' - 7t = 0, 



(6) 



where t = [ti,t2]'^. Here the symmetric 2x2 matrices Sjj, /3jj, and 7jj are 
given in terms of the components of the stiffness matrix C (Destrade, 2001) 
or of the components of the reduced compliance matrix s' (Ting, 2002a) by 
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where X = pc^ and 
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We recall that for two-dimensional deformations of a monoclinic material 
with the plane of symmetry at X3 = involving the coordinates Xi and X2 
only, the relevant non-zero stiffnesses and reduced compliances are related 
through 



(9) 



Finally the boundary conditions are written in a much simpler form than 
when displacement components are involved, as 
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ti(0) = ^2(0) =0, and ti(cx)) = ^2(00) = 0. 



(10) 



3 The secular equation 

3.1 The characteristic polynomial 

Now we seek solutions to the equations of motion 

t(x2) = e'^^'^T, 



in the form 



'11^ 



where '^{p) > 0, to ensure the decay of the wave amplitude away from the 
free surface, and T is a constant vector. So we have by ([6]), 



-anp^ + l^up - 711 l3i2P - 712 



P12P - 7l2 



-a22P - 722 



(12) 



Hence, for nontrivial solutions to exist, p must be the root of a quartic, 
which corresponds to the determinant of the matrix above being equal to 



zero. This quartic, the characteristic polynomial of the equations of motion, 
may be written as 



p'^ + 2u3p^ + U2p'^ + 2uip + Uo 



0, 



(13) 



where the coefficients co^, CO2, coi, and ojq are given by 
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(14) 
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— [4 - ^(s'(l, 2) + s'{2, 6)) + X^s'il, 2, 6)]. 
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Here, the expression s'{ni . . . nk\Tni . . . rrik) represents the determinant of the 
k X k matrix which is a submatrix of the matrix s^- {i,j = 1, ... ,6) and 
whose components correspond to the intersections of the rows ni, . . . , n^ and 
the columns mi, . . . , rrik- Moreover, when ni = mi, . . . ,nk = ruk, the shorter 
expression s'(ni . . . Uk) = s'(ni . . . nk\ni . . . Uk) is used. The quartic (IT3|) was 
obtained by Ting (2002a, 2002b), and by Furs (1997) in terms of invariants 
of the stiffness matrix C. Note that when X = pc? = 0, the quartic of the 
elastostatic case is recovered as (Steeds, 1973, p. 72) 

s'n/ - '^s[y + (2s'i2 + Sqq)p'^ - 2s2qP + s'22 = 0. 



(15) 



Now we use the boundary conditions (TTO|) i 9 at the free surface to estabhsh 
the secular equation. Let pi and p2 be the roots of the characteristic poly- 
nomial (IT3I) with positive imaginary part, and let T*^*") be a vector satisfying 
(I12p when p = Pr{r = 1,2). These vectors are in the form, say. 



T(i) 



a22Pi + 722 

A2P1 - 7l2 



rp(2) 



«22P2 + 722 

A2P2 - 7l2 



(16) 



Then, assuming pi ^ P2, the tractions t defined in (^ are a combination of 

T*^-'^'' and T*^^) for some constants gi and q2. 



t(x2) = gie^'^Pi^^T^^) + gse'^P^^'^T^^). 
Using Ting's (2002a, 2002b) notation, t may be written as 

t(a;2) = B < e^^"* > q. 



(17) 
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where q is the vector [qi, q2]'^, and the matrices B and < e*^^* > are defined 
by 

B = [T(^\ T^% < e*^^* >= diag(e^^Pi^% e'^^P^^^). (19) 

The tractions satisfy the boundary conditions (ITO|) i . that is t(0) = 0, when 

Bq = 0. (20) 

This system has non trivial solutions when detB = 0, that is when the 
following secular equation is satisfied, 

a22^12PlP2 - "22712(^1 + P2) - /3l2722 = 0. (21) 

Now we try to obtain more satisfactory expressions for this equation. 

3.2 The "covariant" secular equation 

First, we decompose pi + p2 and piP2 into their real and imaginary parts as 

P1+P2 = u^ + iu~, piP2 = v^ + iv~. (22) 

It is known that when the roots of the quartic flT5]) are pi, p2, Pi, and P2, 
then ti"*", u~, v~^, and v~ satisfy 

oji = -M+t.+ - u'v~, uo = (f +)^ + (^;")^ (23) 

which leads to the following cubic for v~^, 

{v^f + b2iv+)'' + bi{v+) + bo = 0, (24) 

where 62 = —1^2/2, &i = W1W3 — ujq, and 60 = [1^0(1^2 — 1^3) — '^i]/2- On the 
other hand, the secular equation 02 ip may also be separated into its real and 
imaginary parts, 

"22 A2(^^^) -"22712 (m^) -/3l27l2 =0, «22^12(^^") " «22712(m") = 0. (25) 

At this point, it is important to emphasize that the system of six equations 
(123|) and (125|) for the five unknowns m"*", u~ , f"*", v~ , and X, is consistent. 
This is due a fundamental result of the modern theory of surface waves in 
anisotropic elasticity by Stroh (1962), which states that the complex secular 
equation (ETj) is actually equivalent to a single real equation (for alternative 
proofs see Currie, 1974; or Taylor, 1978). In fact, it can be proved that 
with the appropriate normalization of the Stroh eigenvectors, the real and 



imaginary parts of the secular equation are equal but opposite in sign (Ting, 
1996, p.469). 

Now, v^ must satisfy the cubic (12^ and also the linear equation fl25|) i. 
which using (l23|) i. is 



v+) + ao = 



where a ^^ ^[^6 " ^^U 6|1, 2)] 4, -.-(1,2) 

Working along similar lines but with the displacements components ([T]) rather 
than with the traction components ([5]), Furs (1997) showed that v^ was 
simultaneously the root of a cubic (as in (jMD) and of a quadratic (in contrast 
to 026p ). By writing the resultant of those two polynomials, he obtained the 
secular equation in "covariant form" (his wording) as corresponding to the 
nullity of a 5 X 5 determinant. The same approach applied here to flMl) and 
fl26|) . would yield the "covariant" secular equation as corresponding to the 
nullity of a 4 X 4 determinant. However, v^ is easily deduced from fl2B]) as 
f "*" = — ao and substituted into ^IM to yield the "covariant" secular equation, 

al - b2al + hao - &o = 0. (27) 

Although not stated as such, Furs's "covariant" secular equation is a poly- 
nomial of degree 6 in X = pc^. Equation (127|) is a polynomial of degree 9 in 
X = pc'^, but it may be factorized into the product of two polynomials, one 
of degree 5, which corresponds to spurious roots, and one of degree 4, which 
is the quartic secular equation. This very equation is obtained in a different 
and more direct manner in the next subsection as Eq. (l31l) . 

3.3 The quartic secular equation 

In order to obtain the quartic secular equation directly, we recall that the 
vectors T^^\ T^^\ in fITBl) were computed using the second line of flT^ . When 
the first hne is used, these vectors are 



T(i) 



l3l2Pl^- 7l2 

ttiiP? - AlPi + 7ii 



rp(2) 



l3l2P2^- 7l2 

"iipi - A1P2 + 7ii 



(2^ 



By the same steps that lead from (IT7j) to fl2T|) . we obtain the following alter- 
native form of the complex secular equation, 

anA2PiP2 - "11712 (Pi + P2) + Ai7i2 - /3i27ii = 0. (29) 



By simple comparison of fl^Tl) and fl2^ . the complex terms involving piP2 
and pi + p2 may be eliminated, and the secular equation becomes real, 

aiiA2722 + a22Ai7i2 - a22/^i27ii = 0. (30) 

This equation is a quartic in X = pc^. It was established by this author 
(Destrade, 2001) in terms of the stiffnesses and by Ting (2002a) in terms 
of the compliances (see also Currie (1979) for a less explicit expression, ob- 
tained by writing the equations of motion for the displacements instead of 
the traction components). This polynomial of degree 4 in X = pc^ is also the 
one obtained in the previous subsection by factorization of the "covariant" 
equation fl27|) . Note that it may be obtained directly by using an adequate 
combination of the vectors fITBl) and (1251) for the colomns of the matrix B. It 
is written explicitly with the coefficients in terms of the reduced compliances 
as 

diX^ + d^X^ + ^2^2 + t/iX - 1 = 0, (31) 

where 

_ / / ' _|_ '2 ' I '3 _ / '2 I '2 / \ 
'^12'^11'^22 + ■^le-^ll + -^12 '5ll'5l2 + ■5ll'522j5 

"3 = ~^"5ll"522 ~^ '^11 ~ '^ll'^12 ~r 3SggS]^]^ + S]^2'5ll 

— Z,SiqSii + ^'5l6'5ll'526 ~ ^'^66'^11'^12 + '^22'^11'^12 ^ '^16'^125 
02 = — 3S]^^Sgg + S^Q — 2S]^gS2g — OS^^ + S^yS22 + 2s^-|^s^2 + ■^m^l2i 
"1 = iSii — Si2 + Sgg. 

(32) 



Values of the Rayleigh wave speed cr are given for the 12 crystals of Table 1 
in (Destrade, 2001). They correspond to the square root of the least positive 
root of dSI]). 

In the isotropic case, the reduced compliances take the following values, 

and s'lQ = S2g = 0, where cp and cs are the speeds of the longitudinal and 
transverse bulk waves, respectively. Then the quartic f l3T|) factorizes into the 
product of a polynomial of degree one and of the cubic found by Rayleigh 

(1885), 

(S)'-Kl)^-(--f)S-K-i)-- <-) 



4 Applications 

Now we apply the results of the previous section to other settings, namely, 
to the case of a surface wave propagating in any direction in the plane of 
symmetry X3 = for monoclinic crystals; then to the case of a surface wave 
propagating along a material axis for rhombic crystals; and finally to the 
case where the half-space is made of an incompressible monoclinic material. 
We also seek exact analytic solutions for the speed of Rayleigh waves. 

4.1 Rotation in the plane of symmetry and explicit 
wave speeds 

Monoclinic materials with the plane of symmetry at 0:3 = have a stiffness 
matrix C which is 'structurally invariant' that is, a matrix whose components 
which are zero remain zero after a rotation of the coordinate axes around the 
xs-axis (Bond, 1943). Ting (2000) proved recently that the submatrix of the 
reduced compliance matrix appearing in (Q is also structurally invariant. 
Here, we exploit this property in order to derive the secular equation for sur- 
face waves propagating in any direction in the plane of symmetry. Previous 
efforts covering this topic include those of Chadwick and Wilson (1992) and 
of Furs (1997). 

First, we consider a surface wave propagating on the plane Xg = and 
polarized in the x*-direction, where the coordinate system x* is obtained 
from the material axes coordinate system Xi through a rotation about the 
Xs-axis by an arbitrary angle 6, say. Hence, 



X. 
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0" 




Xi 


—n 



m 





1 




X2 



m 



cos 9, 



n 



sin^. 



(35) 



Next, following Ting (2000), we infer that all the results from the previous 
section are directly applicable to this wave, as long as the reduced compli- 
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ances s'^j, {i,j = 1,2,6) are replaced by the following 'starred' quantities, 

/ / \* / 4 I /o ' I ' \ 2 2 I / 4 I o/ ' 2 I / 2\ 

[s-^^i) = Sj^^m + (/Si2 + -^ee)"^ ^ + ■^22'^ + ^(Sj^gm + 520^ jmn, 

(S22)* = ■§22^'^ + (2s'i2 + -^ee)"^^"-^ + -^ii^^ ~ 2(s2Qm^ + s[Qn^)mn, 

i^'uT = ^12 + (^11 + ^2 - 2si2 - SQQ)m'^n'^ - (s'^g - ^e)!"^^ " n^)mn, (36) 
/ / \* / 4 / 4 0/ ' ' \ 2 2 

— [2s[-j^m'^ — 2s22^^ — (2512 + ■^ee)!"^^ ~ n^)]mn. 



(^e)* = ^26^^ - ^16^^ + 3(si6 - s'^^]m"n 



I \* I 4 /4io/' /\22 

2 r,„/ „2 /o„' I „' N^™2 2\ 



+ [2s22^ — 2s'^]^n — (2s']^2 + ■566)("^ — 72 )]mn, 
{s'mT = s'm + 4(sn + S22 - 2s'i2 - ^e)"^^^^ - 4(si6 - 46)("^^ - n'^)mn. 
In particular, the secular equation for the surface wave is the starred version 
of dSH), that is dlX^ + d^X'^ + d^X"^ + d\X -1 = 0. 

Because the coefficients d\, d^, ^2, and rf* in this quartic are functions 
of 6, it might be possible that for certain angles, the quartic turns into 
a biquadratic, for which the real root X may be found explicitly. Some 
work has been devoted to the search of explicit expressions for the speed 
of elastic surface waves. For instance. Lamb (1904) noted that Rayleigh's 
cubic equation 1^^ factorizes into the product of a polynomial of degree one 
and of a quadratic in c^/c| when Poisson's ratio is |, that is when cj:, = 
3c| or equivalently, when the two Lame constants are equal; in that case, 
c^ = 2c|(l — l/\/3). In the general case, closed-form expressions are rather 
cumbersome for the relevant root of the cubic in isotropic (Nkemzi, 1997) or 
orthorhombic (Romeo, 2001) half-spaces and only approximate expressions 
are sought (Royer and Dieulesaint 1984, Mozhaev 1991). However, Mozhaev 
(1995) showed that for the special orthorhombic materials such that C12 = 
Cge (or with an equivalent relationship for different choices of the material 
axes), the squared Rayleigh wave speed could be obtained as the root of a 
quadratic. Similarly, Ting (2002b) showed that the quartic (1311) simplifies to 
the product of a squared polynomial of degree one and of quadratic in X for 
the special monoclinic materials with the plane of symmetry at X3 = such 
that s[q — 2s2g = s'^2 = 0- Now the starred version of the quartic (^T]) may 
be rewritten in canonical form, 

Y^ + aY'' + bY + e = 0, Y = X + ^. (37) 

Ad*, 



where 



a=[Mldl-?>{dlf]/Wl)% 

h = [(rf*)=^ - Ad;dldl + 8dl{d*^f]/[8{dl)% (38) 

e = [16d;{d;fdl - 3{diy - 256{dlf - 64dldl{dl)^]/ [256 {dl)^]. 
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Clearly, it becomes a biquadratic if 6 = for a certain angle 9 = a say. Then, 
solving the equation 6 = at = a for rfg, for instance, yields the following 
biquadratic. 



Y 

whose explicit relevant root X = pcj^ = Y — d^/{Ad2) is 

d. 



^-~4^\ 



(d-,.2 rff .df.2 1 



Twice, the resolution allowed for a plus or a minus sign. One time, the plus 
sign was selected because X = pc\ must be positive; the other time, the 
minus sign was selected by continuity with the known result (Lamb, 1904) 
in the special isotropic case where Poisson's ratio is |. 

Of course, the existence of an angle 6 = a such that 6 = is not guaran- 
teed. However, numerical simulations show that at least for diallage, gypsum, 
and tin fluoride, two angles a may indeed be found such that a surface wave 
propagating in the x"-direction with attenuation in the a;2-direction, where 
(x^jXg) are obtained from the material axes (xi,X2) by a rotation about 
the xs-axis of the angle a, has a velocity cr = yX/p which is given ex- 
plicitly by (140|) . For diallage, the angles are a = 80.24° and 87.64° with 
corresponding Rayleigh speeds cr = 3960 and 3952 m/s, respectively; for 
gypsum, the angles are a = 18.82° and 65.00° with corresponding Rayleigh 
speeds cr = 2895 and 2946 m/s, respectively; and for tin fluoride, the angles 
are a = 4.01° and 31.21° with corresponding Rayleigh speeds cr = 1324 and 
1351 m/s, respectively. 

4.2 Orthorhombic crystals 

For orthorhombic crystals, s[q = Sgg = 0, so that /^n = 712 = and the 
results obtained for monoclinic materials are greatly simplified. In particular, 
the characteristic polynomial flT3l) is now a biquadratic in p, 



p*-Sp^ + P = 0, (41) 

where the real scalars S and P are given by 

P = ^(l-Xs'J[s',,-Xis[,s',, - s[l)]. (42) 
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Note that, depending upon the sign of S^ — AP and of S, the roots pi and p2 
of the biquadratic with positive imaginary parts are either purely imaginary: 



pi = i\ U-S + V^2_4p)/2, p2 = i\ K-S - V^2_4p)/2 (when S^-AP > 



0, S* < 0), or of the form: pi = a+ib, p2 = —a+ib, where a = \/{S + 2a/P)/4 



and b = y {—S + 2vP)/4 (when S"^ — AP < 0). In any case, we have 

PiP2 = -Vp. (43) 

It is now easy to see that the secular equation (1211) is simplified to (using 

a22VP + 722 = 0. (44) 

Explicitly, the secular equation f l4^ is written as 



(1 - X^)Vl-X46- X^s',,[s',, - X{s',,s',, - s'll )] = 0. (45) 

The secular equation for surface waves in orthorhombic crystals was first 
established by Skelvo (1948) in terms of the elastic stiffnesses. Ting (2002b) 
found the cubic secular equation in terms of the elastic reduced compliances, 
an equation which may be deduced from f H5|) by rationalization; however, the 
squaring process introduces spurious roots, while the exact secular equation 
fH^ has a unique root (Romeo, 2001). 

Note that after rotation about the xs-axis, the secular equation is again 
dlX^ + dlX^ + dlX"^ + d\X -1 = 0, where the d*'s are given by ([32]) and ([MD 
with s'^g = S2e = 0. Also, special directions in which the secular equation 
is a biquadratic might also be found for orthorhombic crystals, following a 
procedure similar to the one exposed in the previous subsection. 

4.3 Incompressible monoclinic materials 

According to Klintworth and Stronge (1990), "many anisotropic compos- 
ite materials appear relatively incompressible because their bulk modulus is 
large compared with their shear moduli. In particular, low-density cellular 
materials are highly compliant in shear because the flexural rigidity of the 
cell walls is small." Nair and Sotiropoulos (1997, 1999) also studied incom- 
pressible anisotropic materials; in particular they considered surface waves 
in monoclinic materials with the plane of symmetry at X3 = 0, but did not 
obtain the secular equation explicitly. 

Recently, Destrade, Martin, and Ting (2002) proved that in linear ani- 
sotropic elasticity, the constraint of incompressibility implied that certain 
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relationships must be satisfied for some compliances. In the present context, 
the following relationships must hold, 

■^11 ' '^12 ~ '^12 ~r -^22 ~ '^16 "•" "^26 ~ ^' \^^) 

and they greatly simplify the quartic secular equation fl?Tl) to 

+ [34 - As[,is[, + s',,)]X' + {As[, + 4)X -1 = 0. (47) 

This equation was obtained by Destrade et al. (2002) in a less explicit man- 
ner. 

Nair and Sotiropoulos (1999) introduced the constants a, (3, and 7 defined 

by 

« = ^^^^, /5 = f^-l' 7 = -4^' (48) 

When these equations are solved for s'^, s[q, and Sgg, the secular equation 
(H7|) may be written as a quartic in x = X/a = pc^/a, 

2x^ - 5(4/3 + 4 - 7^)^^ + (16/3 + 20 - 87^) (4/3 + 4 - 7^)0;^ 

- 4(/3 + 2)(4/3 + 4 - -f^x + (4/3 + 4 - 7^)3 = 0. (49) 

Now a numerical example is given, as the surface wave speed is computed in 
the case (Nair and Sotiropoulos, 1999) where /3 = 0.3 and 6 = 0.1. Then, 
the secular equation (149|) is the quartic 

2.0000x^ - 25.950x^ + 128.56x2 - 247.81x + 139.80 = 0, (50) 

for which the least real root is x = 0.94671, in agreement with the likely limit 
of the 'iterative solution' given by Nair and Sotiropoulos (1999). 
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